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I. INTRODUCTION 



While most working physicists pay homage to the Copenhagen interpretation of the 
Heisenberg-Schrodinger quantum mechanics (QM), many others seek a more causal re- 
interpretation. One ambitious effort in this direction has its origin in the works of David 
Bohm . Bohm employed a formalism for computing the paths for quantum mechanical 



particles close 



y analogous to the method of Hamilton and Jacobi. Nevertheless the Bohm 



approach 0, I^O], also known as Stochastic Quantum Mechanics (SQM), has been thought 
to reproduce in all instances the same probability distributions as does QM. 

The research concerning SQM involves a considerable number of authors dealing with 
various aspects, even if not always within the terms as originally proposed by Bohm 
Studies of extension of the Bohm approach to the relativistic case are also available [?], 

Notwithstanding the relatively difficult nature of SQM computations, it is a 



generally accepted belief that, where a comparison is possible, SQM and QM would give the 
same results. This is indeed the case for the average values of observables at a fixed time. 

However, in a previous work [l2|, examples were reported in which SQM produces results 
different from QM. The examples involved the two-time correlation functions of the electric 
dipole moment components in the hydrogen atom. Within SQM, an explicit numerical 
calculation was performed [12] yielding the related spectral function (Fourier transform of 
the two-time correlation function) for the hydrogen atom in the excited state \nlm) = 1 21 1) . 
It was found to be quite different from the corresponding quantity in QM. 

Our purpose is to study in more detail the comparison between the frequency spectral 
functions in the two theories. The general definition of two time correlation functions (for a 
generic system described by a time independent Hamiltonian H) is 

$ AB (t,t')= 1 -(A(t)B(t') + B(t')A(t)}. (1) 

The related noise spectral functions are defined as the Fourier transform 

VLB/ \ I Jujt^AB/^ al 



S AB (oo) = / e luJt $ AB (t) — . (2) 

J-oo 2tt 

Sum rules for the functions S (u) are investigated in terms of their fc t/l -order moments 

r+oo 

T (*)= / u k S AB {u)du. (3) 



Taking up (for definiteness) the case of the electric dipole moment p(t) = er(t) of the 
hydrogen atom, the following issues are addressed: (i) An explicit analytic form is derived (in 
terms of modified Bessel functions) within SQM for the noise spectral function of the dipole 
moments for the first few excited states, explicitly \nlm) = 1 21 1) , 1 322) , 1 321) and |311). 
General formulas are given for special combinations of quantum numbers, e.g. \n, l,m) = 
\n, n — 1, n — 1) and \n, n — 1, n — 2). (ii) A general proof is given concerning the asymp- 
totic behavior of the spectral functions of the dipole moment fluctuations. As uj — > oo, 
the spectral functions vanish with a power law for both QM and SQM, but with different 
exponents, (iii) The moments of the spectral functions are investigated both for QM and 
SQM, showing explicitly the second order moment differences. 

In Sec. II a brief review of the Stochastic Quantum Mechanics is provided and in Sec. Ill 
the general notion of the two-time correlation function is defined. The related spectral 
function is also introduced. In Sec. IV, an explicit calculation is reported for the hydrogen 
atom two-time correlation function of the electric dipole moment. The spectral function 
asymptotic behavior for large frequency is discussed. Sec. V contains a discussion of moment 
sum rules, and in the concluding Sec. VI the differences between QM and SQM are further 
explored. 



II. THE STOCHASTIC QUANTUM MECHANICS AND PARTICLE TRAJEC- 
TORIES 



In SQM the wave function ip(r,t) entering into the Schrodinger equation, 

H - h i v2+v - (4) 

is conveniently written in the form ip = Re lS / n . R and S are real functions. SQM then 
provides a causal interpretation for the two resulting coupled differential equations thus 
obtained; i.e. 

OS |VS| 2 „ ffi 2 \ V 2 R n 

m + -^r + v -\T»)— = > (5a) 

<w +v .r*™) = o. ( 5b, 



dt \ n 

The first of the above two equations is of the Hamilton- Jacobi form. This is generalized by 
the presence of a new term which takes into account the quantum effects via the quantum 



potential contribution 

, fh 2 \ V 2 R(r,t) 

S is the SQM version of the Hamilton principal function. The equations of motion are 
computed from 

r(t)=v(r(t),t), (7) 



v (r, t) = -VS(r, t) = -Imitp- 1 Vip) (8) 
ji p 

here v(r,t) is the velocity of the particle that passes through the point r at time t. 

Eq. (j5b|) highlights the statistical character of the theory and is interpreted as a conti- 
nuity equation. As with the classical statistical description, one introduces quantities which 
account for the particle properties of an ensemble of identical systems (same Hamiltonian, 
same quantum state, etc.) with trajectories. If the initial distribution in configuration space 
p(r,t ) is assumed to be given by p(r,t ) = R 2 (r,t 0) = \^(r, t )\ 2 then the distribution 
p(v,t) satisfies the continuity equation provided that p(r,t) = \ip(r,t)\ 2 at all times. This 
expresses the time invariance of the configuration-space measure u p(r,t) d 3 r" . The proba- 
bility distribution in configuration space given by p = \ip\ 2 is called the quantum equilibrium 
distribution. A system is then said to be in quantum equilibrium when its configurations 
are randomly distributed according to the quantum equilibrium distribution |4j. This is 
the so-called quantum equilibrium hypotesis (QEH): if a system is described by the wave 
function ip then its configurations are distributed according to p = \ip\ 2 - 

The fact that R 2 (r,t) = \ip(r, t) \ 2 is the probability density that the particle is at r at time 
t holds true in SQM, assures that one finds the same results as in QM. The probability density 
|^(r,0)| 2 gives information on the initial conditions necessary for the quantum Hamilton- 
Jacobi theory to be applied, thereby allowing the determination of particle trajectories 
through Eqs.(J7J). 

In SQM, the probability density at time t is related to particle trajectories. Closely 
analogous to classical statistical mechanics, starting from the initial distribution p(r ,t ) = 
(^(^o, ^o)| 2 , one has: 

p(r,t) = J d 3 r 5[r -r(t,r )}p(r 0l 0) = mr,t)\ 2 . (9) 

Averaging the initial position r*o then yields the same average values of the corresponding 
operators in quantum mechanics. 
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Consider the Hermitian operator A = A(f, p). In the r representation, in a state i[)(r, t) = 
(r\ijj(t)), one has 



(A) t = ip*{r,t)A(r,-ihV)?p{r,t)d 3 r. (10) 



That A is Hermitian, allows the definition of a local expectation value which, when integrated 
over all space, yields the average value {A). One defines 




such that 



(A) t = / A(r,t)R 2 (r,t)d 3 r. (12) 



III. TWO-TIME CORRELATION FUNCTIONS 

In the previous section it was shown how the two theories (SQM and QM) are completely 
equivalent if one considers average values of operators at one fixed time. As already an- 
ticipated above, in order to distinguish the two theories, one needs to consider two-time 
dependent quantities such as the correlation functions defined in Eq.(pQ). The averaging 
procedure is specified in what follows for both QM and SQM. 

A. Stochastic Quantum Mechanics 

As discussed in the previous section, the Bohm approach is able to reproduce the density 
distribution at time t from the initial probability distribution p(r ,t ) = \tp(ro,t )\ 2 at time 
to while still employing the notion of trajectories (see Eq. (jHJ)). For the time evolution of a 
general quantity A, the quantum equilibrium hypothesis (QEH) gives (see Eqs.iPand fT^j) : 



The average value of two operators depending on different times is a generalization of 
Eq.(JT3j). Given any two dynamic variables, A(r,t) and B(r,t), the average value of their 
product, weighed over the initial condition by the weight function \ip(r , 0)| 2 = p(r ), is just 
the SQM two-time correlation function[13]; It is 




(13) 




(14) 



where B(r,t') = B(r(t', r ), 0) and similarly A(r, t) = A(r(t, r ), 0). 
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B. Quantum Theory 



To simplify the notation, consider the case of a time independent Hamiltonian with a 
discrete spectrum of eigenvalues 

H\N) = E N \N). (15) 

In QM the two-time correlation function for a given state \N) is 

K B (t - = ~ (iV| A(t)B(t') + |iV) , (16) 
where (in the Heisenberg representation) 

A(t) = e iHt/n^ e -im/n_ ( 17 ) 

Consider the special case A = B. Then Eq. (JTBJ) reduces to 

H A (t -t') = \ (Ki(t - + - 0) , (is) 

where 

$^(t ~t') = ^J2\ ( M l ^ 1^) 2 e™^^*') (19) 

and Hujmn = (Em — En). The Hamiltonian being time independent is reflected by the fact 
that the $ N ' B depend just on (t — if). 

The related noise spectral functions, defined as in Eq.(j2J) in QM are easily found to be 
given by 

S N (u) = S^(u)+Sjf(u), (20) 

with 

A/ 

When |7V) is the ground state only S N contributes for uj > 0. In general, Sn(lj) is an even 
function of uj. Furthermore, if the Hamiltonian has a spectrum with both a discrete and a 



continuous part, then the sums in the above Eqs. (j2()J21j) clearly split into a corresponding 
sum over the discrete part of the spectrum plus an integral over the continuous one. It will 
prove useful in deriving the asymptotic behavior of the QM noise spectral function Sn{uj) 
(in the next section) to relate it to the Fourier transform of the average of the time-ordered 
product; i.e. 

? r +OQ r * - i 

a N {u) = -\ e mt {N\T A(t)A(0) \N) dt, (22) 



3C 



where "T" denotes operator time ordering. The quantities Sn(u) and axiyj) are related by 



Sn(w) = 7T Im «jv(w). (23) 
Zix 



IV. THE ELECTRIC DIPOLE MOMENT AND THE HYDROGEN ATOM 
A. Stochastic Quantum Mechanics 

Consider the electric dipole moment along the x axis d x (t) = ex(t). In the Bohm descrip- 
tion, the two-time correlation function of the electric dipole moment for a generic excited 
eigenstate of the hydrogen atom with quantum numbers (n, I, m) is as given in Eq. (|14|): i.e. 

<WM') = (d x (t)d x (t')) nlm 

= J d x (r(t,r ))d x (r(t',r ))p nlm (r )d 3 r , (24) 

where d x (r(t, r*o)) = er sin 9q cos 4>(t). Using the hydrogen atom bound state wave functions 
and the solution to the equations of motion in the Bohm description, given in appendix A, 
it is straightforward to derive the following properties of the related noise spectral functions 

S nt i j0 {u>) oc 5(u). (25) 

Thus, without loss of generality, m > and u > can be assumed in the following discussion. 
For general quantum numbers, one can derive: 



(26) 



where 



= (h/pa 2 



0J> 



2_ jmujo 



I1V U) 



2 



(27) 
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n\2l + m)!((2m - l)!!) 2 (n -/-!)! 



2nm(l + m)\[{n + Z)!] 3 



and Cj^ 1 ^^) are the ultra-spherical Gegenbauer polynomials |14| which satisfy the follow- 



ialsQ 



(28) 



ing relations: 



^(m+1/2) 



J l—m 



a 



(m+1/2) 




(1) 

(0 



(I + m) 



(2w)!(Z - m)\ 



(29) 



1-3- ••(2m- 1) 

where Pi(£) is the Legendre polynomial. 

A few special cases have been explicitly computed and are here reported using S n i m (u) 
to denote the noise spectral function of Eq. (j26|) in units of (e 2 al)/u>o : 



cSQM/ s 



O n i m (UJ) 



n = 2 
n = 3 

5322 (<^) = 
<S32lM = 
S3I1M = 

n (i = m = n — 1) 



CJ 



W") = ^ (^) 4 * 2 ,l*l(*2, 



(30) 



1 



w 



2187/ Vcj 



1 



UJ () 



3888/ V ^ 
1 \ (Uq 



^3,2-^1 (^3,2) 

2 3 ,1 [2^1^3,l) + *3,1 ^0(23,1)] 



243 J \w 



— Z 3,l 



5 z 2 \ 7 



8 16 



Sn,n-l,n-l{w) — C n ( J Z n) n-\K\{z n)Tl _\) 



16 



(31) 



2;i 



[(2^-3)!!] 2 
2n[{2n - 2)!p 



(32) 



n (I = n — 1, m = n — 2) 



to 



n+l 



Z n,n-2^2{z n ,n-2) 



Cn,n-l,n-2[(2n-l)!(2n-3)!!] 



2 4(n - 2) 
128n 2 



(33) 
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~i — i — i — i — i — i — i — i — r 



1 — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



In the above expressions, K n (z) are the n th order modified Bessel functions. It should be 
noted that the function SWity) agrees exactly with the numerical computation previously 
reported 3|. It has been verified that the particular results of Eqs.(|HU|) and (j3Tj) are 
consistent with general formulas given in Eqs. (|32|) and (|33j) . In Fig.l we show the plots of 
some of the above explicit examples. 

10 4 
10 3 
10 2 
10 
1 

KT 1 




j i i i i i i i i_ 



0.025 0.05 



0.075 0.1 



0.125 



FIG. 1: Shown are four plots of typical spectral functions S n i m (oj) = ujqS^^ (uj) / {eao) 2 for the 
electric dipole moment of the Hydrogen atom computed employing the Stochastic Quantum Me- 
chanics. 



Asymptotic behavior. 

From Eq. (|2rj|) it is possible to derive for general quantum numbers the asymptotic behavior 
at large frequencies of the SQM noise spectral function; it is 



cSQM/ 

°nlm \ 



C 



nlm 



eV 



o \ ("o\ 



3+m 



as UJ — > OO, 



(34) 



where 



C 



nlm 



128 



2 f Am\ 



3+m 



-(21+1) 
J n+l 



(P) 



dp. 



(35) 
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Thus, noise spectral functions vanish as uj — > oo with a power law. The exponent is 
related to the state's quantum numbers by (3 + m). The explicit cases considered above can 
easily be shown to agree with Eq. (J34|) . One may use the expansion of the modified Bessel 
functions for small values of the argument |l4j, i.e. as z — > 



KJz) 



(l/2)r( I /)(V2)-" if p^O, 
-In(z) ifz/ = 0. 



One then infers the large u — > oo behavior of the spectral functions; i.e. 



qSQM/ 



(1/12* 



for state |211) 



1/2187) f — ) for state |322) 



a; / 



(1/1944) (— ) for state |321) , 
(5/1944) (— ) for state |311) . 



(36) 



(37) 



B. Quantum Mechanics 

As anticipated at the end of Sec. Ill, in order to derive the asymptotic behavior of the QM 
noise spectral function it proves useful to work with the "time ordered propagator" defined 
in Eq.(j22J) and use Eq.(j23J) to find Sn(lj). The function a^{oj) can easily be connected with 
the (retarded) Green's function of the Schrodinger equation; i.e. 



• poo 

G ret (r 2 , r i; u) = -- J (r 2 | e^"™* |n) dt. 



(38) 



The time ordered response function <x£ (oo) is related to the Green's function G ret (r 2 , rj.; to) 
via 



-a K N >{uj) 



J d? ri j d 3 r 2 (N\ A(0) \r 2 ) G ret (r 2 , n; ^ + (n| B(0) |JV) + 



d 3 n / d 3 r 2 (N\B(0)\r 2 ) G ret r a ,n; 



IV 



w (n|A(0) |AT> 



(39) 



The retarded Green's function associated to the Hamiltonian if, as defined in Eq. (j38j) . 
satisfies the differential equation 



(hu - H (r2) )G ret (r 2 , n; w) = 5(r 2 - ri), 



(40) 
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where H, 



(r 2 ) 



— (h 2 /2fj,)Vf s — (Ze 2 /r 2 ) denotes a one electron Coulomb system. The 



differential equation o 
malization) is 



(r 2 ) 

the non-relativistic Coulomb Green's function (in the standard nor- 



V (r-2) 



2ku 



k 2 } G(r* 2 , ri; uj) = 5{r 2 - r% 



(41) 



where k = \J (2fiuj/h) and v = AnTi 2 /(Z^e 2 k). 

The Green's functions appearing in Eqs.(jBHJ) and (j41)l are related by a normalization 
constant G ret (r 2 , r%] uj) = (2fi/h 2 )G(r 2 ,ri,uj). A closed expression of the Coulomb Green's 



function in terms of Whittaker functions 

r(i 



G(r 2 ,ri;w) 



IV) 



det 



14j has been given by L. Hostler [l5| as 



(42) 



l7r l r - r > I \W lu ., 1/2 (-tk ai ) M lU]1/2 {-ika 2 ) / 
where the dots over the Whittaker functions denote differentiation with respect to their 



arguments and 



«2 = r 2 + ri + |r 2 - r*i| 
«i = r 2 + r x - |r 2 - rJ 



(43) 



Asymptotic behavior. 

The aim of this section is to derive the asymptotic form of the noise spectral functions 
as u) — > oo. Using Eq. ()39)l requires the Coulomb Green's function in the regime ui — > oo 
(respectively uj — > — oo). In this limit, k — * oo (respectively A; — » zoo). Also z/ — >• so that 
Eq.(|4*T|) reduces to the differential equation of the free particle Green's function Go(r 2 , r\\ uj). 
Indeed, from the exact solution in Eq. (j4~2")) . it is possible to show explicitly that as \uj\ — > oo 

gifc|r2-ri| 



G (r 2 , n; uj) -> G (r 2 , r x ; w) 



(44) 



4vr|r 2 — t*i I 

This proves that the exact Coulomb Green's function has an oscillatory behavior at large 
positive frequencies and an exponentially damped behavior at large negative frequencies: 



— An\r 2 — t*i | G 
— 47r|r 2 — 7*1 1 G 



,l j"2 —ri 



yj 2flU)/H 

l\y/2fi\w\/H 



as u; 



as uj 



+oo, 
-oo. 



(45) 



Inserting this result into the central Eq. (|H9j) . one finds that only the first term survives; i.e. 
for uj — >• oo 



(h 2 /2rr^ B \uj 



„ik\r2— n| 

rf 3 r! / d 3 r 2 (A| A(0) |r 2 ) 1 f (r 2 | B(0) |A) 

r* 2 -7*i 



(46) 
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where it is to be recalled that k = \J (2/jIJJWj. 

The above considerations can be readily applied to the case of the electric dipole moment 
in the hydrogen atom with A(0) = B(0) = ex. As u — > oo 



UP 2 r pik^-r^ 



2-kTi J \r 2 - ri| 

Using the above in Eq. (J23|) . the asymptotic expression for the electric dipole moment spectral 
function is derived; i.e. 

S N X) (") - A 2 \ Im I ^ I d 3 r 2 r N (r2)x 2 f r2 n , x 1 ^(r 1 ). (48) 

With | AT) = \nlm), the hydrogen atom wave functions are written as ip n im{f) = 
Xni( r )Yi m (9 , 4>) . One may also employ the expansion 

f 1 ^ = (Am^f^Mkr^ikr^ V Y^faMYUhM, (49) 

where r < = min(ri,r 2 ) and r> = max(ri,r2). The noise spectral function then has the 
asymptotic limit 

e 2 f t r°° poo . 

S him( u ) 2 Jm 1 i/c $Z / / rfri *" 2 ( r 2ri) 3 X„«(r 2 )jz'(A;r < )4 1) (^ > )x„/(n) k 

(50) 



with the constants C /m defined as 



t(0 ui . ju 

m'-Z' J 



y,,* m ,(0, </>) sin cos (j)Y lm (9, (f>)dQ 



When taking the imaginary part in Eq. (|5()|) . only the function h\}\kr) is complex |l4l|: 



2 



(51) 



Im 



ihf{kr)} - j,(fcr) = v /n/(2kr)J l+1/2 (kr). (52) 



Therefore, the noise spectral function as uj —>■ oo reads 



(2 \ oo 



r 5/2 Xm(r)Ji'+i/ 2 (kr)dr 



(53) 



The radial integral in the above expression can be evaluated in the limit of high frequencies 
and is found to vanish as (/ca )~ ( - 4+ ' +1 ^ 2 - ) . Indeed, using the hydrogen wave functions reported 
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in the appendix, one finds 
Iniv — 



^ / r x 2 



r — Xm(r)J u+ i/ 2 {kr)dr 
o V a o/ 

(n-/-l)!^ 1/2 



(54) 



n 2 V [(n + /)!]* 



(fca ) (4+z - 1/2) V™ / io 



( ^ ;;£ , ;:c (3+«-l/2) e -(x/nfeao) x 



^S 1 {2x/nka ) J v +\/7,(x) 



When k — > oo the Laguerre polynomial can be replaced by the constant value that it takes for 
a vanishing argument L^l(0). The remaining integral is tabulated (see Eq.(6.621) in [18.]) 



and one finds 



L nlV 



(n-l-l)\ 



1/2 



n 2 V [(n + Z)!] 3 J (/va ) {4+ '~ 1/2) \n 



L^(0)r(l + l' + 4)x 



P. 



-(/'+1/2) 



3+2-1/2 



(0) 



nkao 



dx ^3+2-1/2 W 



(55) 



where P^(x) are the associated Legendre functions of the first kind. One should note that 
for the sum in Eq. (j53j) only few terms are non- vanishing. This is related to well known 
electric dipole selection rules which apply when calculating the quantities 

i' 

Cf2= W' m'\ sin 8 cos <f)\lm)\ 2 , 



m'=-V 



i.e. 



(l'm'\ sin 9 cos <j) \ lm) ^ only if 



m 



mil 



V = l±\ 



(56) 



Applying the selection rule, I = I ± 1, it turns out that -P 3+ ;_ 1//2 (0) vanishes for any I, 
while its first derivative at zero is always finite (see Eqs.(8.6.1) and (8.6.3) in Q]). This 
completes the proof of the asymptotic behavior of the noise spectral functions in QM. Taking 



into account that ka$ = ^J2uj/ujq, one finally concludes that for uj — > oo, 

-QM, , / < "n \ ,,, fWo\ 4+l+1 / 2 



I ^nlra [ 
U>0 J \w 



where 
C 



nlm n 4 



£ c«r<^i^[LS'(o)] 2 



2'=2±1 



[(n + /)!]3 



r(z + r + 4) 

n(v / 2) 4+/+1 / 2 



d 



p-('' + l/ 2 )/ n N 

dx w - 1 ' 2 1 j 



(57) 



(58) 



When comparing Eq. (|57|) with Eq. (|34jl a difference in the two descriptions, SQM and QM 
is made very clear. As uj — > oo, SQM predicts for noise spectral functions in state \nlm) to 
decrease at large frequencies as u~^ 3+m ^ and QM as uj~~^ +1+1 / 2 \ 
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V. SUM RULES: MOMENTS OF THE NOISE SPECTRAL FUNCTIONS 



In highlighting possible differences between the predictions of Quantum Mechanics and 
the Stochastic Quantum Mechanics, it proves useful to study some global properties of the 
spectral function, e.g. sum rules. The zeroth order moment is readily evaluated, 



hoc 



7 (0) = / S A ' A {u)du 

dudt 



eo 
+00 



e lU}t G A ' A {t)- 



2tt 

= G A ' A (0) = (A 2 (0)) . (59) 

The zeroth moment 

7 (o) is 

nothing more than the average value of (A 2 ) at time zero and 
it is the same in QM and SQM. For the special case under consideration this assumes the 



value (x 2 (0)) 



nlm' 



7il = ^[5n 2 + l-3/(/ + l)] 



(l + m + l)(l -m + i) (l+m)(l-m) 



(60) 



(2/ + l)(2/ + 3) (2/ + l)(2/-l) 

On the contrary, the calculation of the second order moment 7^) (again in both theories) 
shows the first discrepancy between QM and SQM. The computation of the second order mo- 
ment can be related to the second derivative of the two-time correlation function calculated 
at r = t — t' = 0. Indeed repeated integration by parts yields 



-00 



u S N {u) = - e - — . (61) 



00 



dr 2 \2n J 



Thus 



Tn = u S N {uj)duj = - — — — . (62) 



dr 2 



T=0 



A. Stochastic Quantum Mechanics 

In the Stochastic Quantum Mechanics from Eqs. ()24j) and (|A6jl one has 

dr 2 \ fi J r 2 sm 2 9 



so that 

(2) e 4 /m 



(64) 
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In QM one obtains 

d 2 $ni m {r = 0) e 2 



(nlm\{x(0),x(0)}\nlm) . (65) 



dr 2 2 

Using the Coulomb Hamiltonian H = (p 2 /2yu) — (e 2 /r) one may deduce the following com- 
mutators: 

x = -[H,x\ = — , 
n n 

i r„ , f e 2 \ x 



so that 



* = s^J = -UJ^ (66) 



d 2 $ni m (T = 0) e 4 x 2 

(nim\ — \nlm) (67) 



and 



X 2 ) 



7(QM)nZm ^ ^ n 2 



1 



dr 2 2/x r 3 

(Z + m + 1)(/ - m + 1) (l + m)(l-m) 



(68) 



(2/ + l)(2/ + 3) (2/+ 1)(2Z — 1)_ 

which is quite different from the expression obtained in the Bohm theory Eq. (j64|) . Clearly the 

differences found in the second order moment imply rigorously different spectral functions. 
B. Semi-classical limit 

It is interesting to see how the two quantities which are different for general quantum 
numbers have the same semi-classical limit for large values of the quantum numbers n, I and 
m. Setting the maximum orbital momentum I = m = n — 1 one finds that 

7 (2) - ( —\ {n ~ 1] (69) 

7( S QM)n,n-l,n-l " [ ^ j ^3 ' 



T(QM)n,n-l,n-l ( o„ ,. L/ n .i/9l ' ( 70 ) 



-A 



2ao/i / n(n + 1/2) 
In the limit of large values of n the two theories agree; i.e. 



lim KtS-^-i} = (J-) (71) 

n->oo \ZaofX J 



for both QM and SQM as expected. 



VI. CONCLUSIONS 



In this work the predictions of the Stochastic and conventinal Quantum Mechanics have 
been compared in some detail with respect to the two-time correlation functions. The 
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example of interest in this work is the dynamic evolution of the electric dipole moment within 
the hydrogen atom. Previous numerical computations [l2j of the SQM spectral function for 
the first excited state |211) have been confirmed. In addition closed expressions have been 
obtained in terms of modified Bessel functions K n (z), for several excited states with special 
combinations of quantum numbers. A derivation is provided of the asymptotic form (as 
uj — > oo) of the noise spectral functions for both QM and SQM. For large frequencies the 
two descriptions provide different power law behavior. For the hydrogen atom bound states 
\nlm) the SQM spectral functions scale oc uj~^+ m ) as opposed to the QM spectral functions 
which scale oc lo~^ +1+1 / 2 \ 

The difference in the noise spectral functions is reflected in different sum rules which are 
obeyed by the spectral functions. Explicit exact evaluations of the second order moment 
of the spectral functions 7^, as in Eq.Q, show that the two descriptions predict indeed 
different values. 

Finally by considering, for example, the interaction of an hydrogen atom with the field 
of an electromagnetic wave, it is possible to relate the noise spectral functions within an 
excited state to a total absorption cross-section by the relation a^ ot {u) = 87r 2 a QED [uj Sn{uj)] 
as shown in detail in the appendix C. 

Let us briefly comment on the different predictions just so pointed out between conven- 
tional quantum mechanics and trajectory based interpretations of stochastic nature (Bhom 
or Nelson). It is certainly worthwhile to mention here a recent work where the author 
shows that trajectory based interpretations of quantum mechanics are incomplete. This 
happens for systems with unbounded Hamiltonians. In particular it is shown that for par- 
ticular systems (providing explicit examples) there exist states of finite energy for which 
the decomposition of the Schrodinger equation into a continuity and modified Hamilton- 
Jacobi equation is impossible. These examples are also shown to be connected to the fact 
that the corresponding state wave functions exhibit fractal properties. The main conclusion 
of ref. Q| is that Quantum Mechanics goes where trajectory interpretations do not follow 
despite their (in principle) duty to do so. So one might wonder whether the results found 
in the present work regarding two-time correlation functions are to be ascribed to such in- 
completeness of stochastic approaches to quantum mechanics. We can just remark that the 
states considered here to evidence differences in the predictions of the two theories are the 
bound states eigensolutions of the Hydrogen atom hamiltonian (\ip n im)) and as such do not 
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have the properties required in |19( to highlight the supposed incompleteness of Stochastic 
QM (i.e. undefined Hip, but with finite average energy). It would certainly be interesting 
to consider the possibility to construct such states for the Hydrogen atom but this deserves 
further investigation, and goes beyond the scope of the present work. Were it possible to 
confirm this connection it would leave little doubt on the authors' minds as to which of the 
two theories would have to be ruled out. 
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APPENDIX A: THE BOHM DESCRIPTION OF THE HYDROGEN ATOM 

A stationary eigenstate of the Schrodinger equation for the hydrogen atom is written as 



^nim(r,t) = X m{r)Y lm {d,<P)e~ iE - t/h 



(Al) 



where 



Xm{r) 





'(n-i-iy: 


(n 2 ) 





1/2 



Y lm ( 



p = (2r/na ) 



N lm P! ml ( cos 6)e im4 



(A2) 



L n +i (p) are the associated Laguerre polynomials and P\ (cos 6) are the associated Legendre 
functions |l4j|. Eq. (jAl|) can then be written as 

VwM) = N lmX ni{r)P} m \cose)e^ nm ^ t \ (A3) 

where Ni m , Xn,i(f) and P ; ' m '(cos#) are real. Thus, the dynamics of the system in the Bohm 
description is provided by the quantum action 



S(r, 9, 0, t) = hmcj) - Et. 



(A4) 
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The Eqs.(J7J) of motion are 



v r = r = - ( ^ ] = 0, 



/i \dr 

v <i> = ( jrr ) = rsin(0)0 

/xrsmd \o<p J 

mh 



(A5) 



fir sin 9 

These can be integrated yielding 

r(t)=r , 9(t)=9 

and 

m = *° + \jfi*?i;)- (A6) 

APPENDIX B: SECOND ORDER MOMENT IN THE STOCHASTIC QUAN- 
TUM MECHANICS 

Let us here consider the calculation of the moments in SQM. According to the definition 
given in Eq.© for the state |211) (and for even n) one has 

2„2 \ poo 



64ujq J Jq \uj J 
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Using the change of variable x = ^uj /uj, one finds 

7&, = -o [ x^ Kl (x)dx. (B2) 

This yields 

7sqm = 12 (ea ) 2 

APPENDIX C: HYDROGEN ATOM INTERACTING WITH THE FIELD OF A 
PLANE ELECTROMAGNETIC WAVE 

Let H denote the Hydrogen atom hamiltonian and suppose to have an atom in one of 
his stationary eigen-states (\N)) at t = to 

!*(*>)> = W) (CI) 
18 



interacting with the electric field of a plane wave of frequency 00: 

V int {t) = E(t) ■ d 
E(t) = E e x cos(kz - ut) (C2) 

the elctric field component of the plane wave assuming that it is traveling in the z direction 
with momentum k = uj/c and d = ex is the dipole operator. 
The full hamiltonian is therefore: 

H = H + V mt (t) (C3) 

The quantity we would like to study is the total transition probability per unit time given 
that the system is initially in the state \N) at time t = t Q . Let us compute first the total 
transition probability py^ arLSltlon ^) a t time t, given the initial condition in Eq. IC11 At time 
t the system will be in the state \^{t)) obtained from the state \N) by application of the 
evolution operator (i.e. solving the Schrodinger equation for Thus the probability Pjv(i) 
that at time t the atom is still in the state |7V) is given by: 

P N (t) = I (N\*{t)) | 2 (C4) 

Conservation of probabily requires that: 

P N (t) + W^ ansit[on \t) = 1 (C5) 

In this problem the interaction has an explicit time dependence and so it is useful to resort 
to the interaction representation: 

= U(t,t )\*(t )) 
U(t,t ) = Te^^'^ 

Kt(t') = e+^'vUty-^' (C6) 

where T stands for time ordering. 

Thus the probability of being in the state \N) at time t is given by: 

P N (t) = \(N\U(t,t )\N)\ 2 (C7) 

In second order perturbation theory the amplitude of remaining in the state |7V) at time t is: 



(y\u{t.to)\x) = i-i f dt'(N\v; nt (t>)\N) - f f dt'dt>\N\T[v; nt (t')v^{t")\ \n) 



f, mtv n 1 on 2 , 

11 J t *n j to j to 



(C8) 
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The dipole interaction we are considering does not contribute at first order since (N\d\N) = 
for any eigenstate of the hamiltonian Hq (hydrogen atom). Thus: 

(N\U(t,t Q )\N) = 1-Ejv 

P N {t) = 1 - 2$(£*) (C9) 



Comparing this last Eq. with Eq. IC5I one finds for the total transition probability 

transition = 2 ^(Eyv) = 21m (lE N ) OV. 

WjT^it) = Im { + J-2 I] I dt'dt"(N\T [VUt')Vl t {t")\ \N) } (CIO) 

Now let us insert the explicit form of the dipole interaction given in Eq. IC2I The fact 
that the wave is assumed to be x-polarized selects the x component of the dipole operator 
e x ■ d = d x = ex and: 

V- nt (t) = eE x(t) cos(kz — ui) 

x(t) = e + * Hot x(<S)e-* Hat (Cll) 

Inserting the above expression of Vi n t(t) in Eq. IC10I we also adopt the so called long wave- 
lenght approximation which consists in neglecting the z— dependence in the interaction. 
This is justified so long as ka <C 1 being a the Bohr radius. Thus one gets: 



^transitions = J m ^ + _ ( e £ Q )2 / / tf, (jy-jy f^')^")] | jV) CO s{ut') COs{ut") 

J to J t 

(C12) 



t rt 



Then note that defining r = t' — t" and t + = t' + t" one has: 

cos(atf') cos(u;t") = - [cos(cjt) + cos(c<ji + )] 
(N\T[x{t')x(t")\\N) = (N\T t [x(t)x(0)]\N) (C13) 

where T T stands for time ordering relative to the r variable. Then define to = —T/2 and 
t = +T/2, change the integration variables according to dt'dt" = (l/2)drdt + to obtain: 

transition ( T ) = j m | + J_ (^o£ ^ + ^ ^ + [ x ( r ) x (Q)] | [cOs(wt) + COs(wt+)] J 

(C14) 

Now the time integration over the variable t + is readily done: 

^ ransition (T) = Im j+i J +T dr (N\T T [x{t)x{0)\ \N) 2T cos(cur) + i sin f ^ | 

(C15) 
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And the total transition probability per unit time y;^ ansitlon (j n the limit of infinite times) is 
extracted: 

Ti/transition/'TA ( • (pip \2 r+oo ^ 

<" Sition =lim N T { > =Im J dr(N\T T [x(r)x(0)]\N) cos(^r)| 

(C16) 

Finally it is easily shown that: 

(N\T T [x(-r)x(0)] \N) = (N\T t au [x(t)x(0)] \N) (C17) 

and hence: 

^transition = j m ^ + ^_ i^f J °° ^ [ x (r) X (0)] |JV) J (C18) 

Thus defining the time ordered "propagator" by: 



+00 



a N (uj) = - I dre +WT (N\T T [x(t)x(0)} \N) (C19) 



one writes the total transition probability per unit time as: 



w 



transition _ ( eE o) 2 /m ( C20 ) 



JN ~ 2h 

This total transition probability when normalized to the flux of incident photons defines 
a total transition cross-section a^ t (uj) (adsorption and possibly stimulated emission, if the 
state \N) is an excited state): 

w 

or : 



/, ,\ „,, transition f f^oi \ 



2 

= &K~^Im[a N (u)] (C22) 

On the other end we have shown, c.f. Eq. [231 that the imaginary part of aj^{uj) is directly 
related to the spectral function (fourier transform) of the two-time correlation functions: 

h r +co 1 

—Im[a N (u)} =S N (u) = / dre + ^ - (N\[x(r)x(0) + x(0)x(r)]\N) (C23) 
^ J-00 2 

We therefore conclude: 

a? ot (u) = 8vr 2 a QED [uS n {cj)] (C24) 
a QED = e V(^ c ) ~ 1/137 being the fine structure constant. 
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